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Abstract 



If light sterile neutrinos are needed to understand the neutrino puzzles, as 



CN . is currently indicated, a major theoretical challenge is to understand why its 



mass is so small. It is a more serious problem than understanding the small 

^ I mass of the familiar neutrinos. We discuss a new way to solve this problem by 

1^ ' identifying the sterile neutrino as gauge neutral fermion propagating in the 

'!;;j- ■ bulk of a higher dimensional theory, with bulk size of order of a millimeter. 

^^ . The smallness of its mass is then a consequence of the size of the extra di- 

Qv^ , mension and does not need the introduction of new symmetries. We present 

^N I a realistic model for neutrino masses and mixings that implements this idea. 

iu: 

(D . I. INTRODUCTION 

Recent experimental data on atmospheric neutrinos from the Super-Kamiokande |ll| col- 
^ ■ laboration has provided conclusive evidence for neutrino oscillations and confirms earlier 

indications of such oscillations from other experiments 0. In addition there is also evidence 
for neutrino oscillations from the solar neutrino deficit observed by Kamiokande, Homestake, 
Gallex, Sage and Super-Kamiokande 0,^, and direct observation of u^ -^ z/g (z/^ — > Ue) in 
the LSND |^ experiment. All these suggest nonvanishing masses for at least two of the three 
standard neutrinos and opens the window to explore physics beyond the standard model. 

It has been realized that to explain all three evidences for neutrino oscillations three 
different mass differences are needed, while only two can be obtained with three neutrinos. 
Solar neutrino data requires ^ 

Am^ ~ 3 X 10"^ - 1.2 X 10"^ eV^, 
sin^ 26^01 ~ 3 X 10~^ - 1.5 x 10"^ (1) 
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for the small mixing angle MSW 0] solution. Atmospheric neutrinos seems to prefer z/^ — z/^- 
oscillations with a maximal mixing |^: 

Amlt^ ~ 4 X 10-^ - 5 X 10"^ e1/^ 
sin^ 29atm ~ .76 - 1. (2) 

Finally, the LSND results along with other constraints from KARMEN 0, Bugey |T^ and 
E776 at BNL W^ suggest 



Amt^^.2-2eV\ 



sin^ 2^eM ~ 2 X 10"^ - 4 x 10"^ (3) 

A simple explanation of these parameters could be obtained if there is a sterile neu- 
trino |]12[ (which is a light neutrino that does not couple to the Standard Model particles). 
The existence of a light sterile neutrino however poses a major theoretical challenge. It is 
more problematic to understand its small mass than to understand the small mass of the 
familiar neutrinos. This is because, the familiar neutrinos are standard model gauge nons- 
inglets; therefore their mass can at most be of order of the weak scale. On the other hand 
the sterile neutrinos being standard model singlets, their apriori mass could be of order of 
the Planck scale- some sixteen orders of magnitude higher. 

Several scenarios involving sterile neutrino have been studied in the literature [[l^ T6[ 



where the general strategy has been to invoke new symmetries of different kinds to explain 
its small mass. It is clearly too early to say which if any of these approaches is the correct 
one. Here we explore a completely different approach to its mass. We consider the possibility 
that the sterile neutrino is a massless Dirac bulk fermion in a higher dimensional brane-bulk 
picture with at least one large extra dimension. The bulk fermion has Kaluza-Klein (KK) 
excitations whose spacing is given by the inverse size of the bulk. If the bulk size is large 
(say a millimeter) as has been speculated recently [jl^, then the sterile neutrino (i.e. the 



bulk neutrino) mass is automatically small without the need for any extra symmetries. In 
recent days, interest in such theories (i.e. theories with large extra dimensions) has been 
heightened from independent experimental considerations since they lead to deviations from 
Newton's inverse square law of gravity at small distances. In fact the maximum bulk size 
believed to be at the boundary of gravity experiments is a millimeter, which corresponds to 
sterile neutrino masses of milli-eV's. This value of masses is precisely what is interesting for 
solving the solar neutrino problem providing an intriguing connection between the neutrino 
experiments and those searching for deviations from Newton's law in the millimeter range. 

In the brane-bulk models with large extra dimensions, the smallness of the KK excitations 
of the bulk fermion mass is related to the size of the extra dimensions as just explained. On 
the other hand the small mass of familiar neutrinos arises from the following fact. There is 
a suppression of the coupling between neutrinos living on the bulk with those fixed on our 
four dimensional world |lT8|-pH] which arises naturally from normalization of the bulk field. 
As far as the oscillations between the brane and bulk neutrinos is concerned, it is dominated 
by the mass of the lowest KK mode of the bulk neutrino. Thus one has the right parameter 
range to give an explanation for the solar neutrino deficit. 

Two distinctive scenarios which involve bulk neutrinos have been outlined in recent 
literature. The first uses only the standard model left handed neutrinos in the brane and 



one or more bulk fermions [18,19| and a second one that has both the left and right handed 



neutrinos in the brane in combination with a bulk neutrino 211. In the second case, one 



could invoke the seesaw mechanism |^ in the brane to understand the smallness of the 
neutrino masses; therefore the higher dimensional physics need not play a role. Of course it 
could happen that implementation of seesaw mechanism is impossible due to group theory 
of the model |2l| in which case the extra dimension as well as the string scale will play a 
role in neutrino physics. 

It must however be emphasized that even though mechanisms to understand small neu- 
trino masses have been outlined, no attempt has been made to construct realistic models 
that explain neutrino observations using the property of higher dimensional physics. It is 
the goal of this paper to make an attempt in this direction and in particular construct a 
model involving a sterile neutrino that explains the present neutrino data. 

We first study the models without right handed neutrinos and point out that even though 
the mass of the sterile neutrino is naturally small in these models, several unnatural assump- 
tions are needed if we are to obtain a desirable mixing and mass pattern for understanding 
observations. We then study models with both left and right handed neutrinos in the brane 
with seesaw mechanism that explains the small masses for two of the known neutrinos and 
include a bulk neutrino with one extra dimension being in the submillimeter range to play 
the role of the sterile neutrino. Since the masses of many of the low KK excitations are 
then in the range of interest for solar neutrino oscillations, the key question is whether the 
mixings pattern is such as to be of interest in solving the neutrino puzzles without at the 
same time contradicting known observations. In particular, we explore whether only the 
lowest mode of this bulk neutrino can play the role of the sterile neutrino. The advantage of 
this approach over other models for the sterile neutrino is that smallness of its mass is now 
a geometrical rather than a symmetry effect. The nontriviality of the problem arises from 
the fact that the bulk neutrino of course has infinite number of excitations and one has to 
tackle the mathematical problem of extracting the physical masses and mixings of neutrinos 
from this complex situation and study whether one has a solution to the neutrino puzzles. 

Our basic results in the second case are the following: the lightest eigenstate which is 
predominantly the electron neutrino is massless and it mixes with the first excited KK mode 
of the bulk neutrino, which has mass of order R~^. We use this to solve the solar neutrino 
deficit via the small angle MSW mechanism. This requires the size of the extra dimension 
to be in the micro to millimeter range and a string scale which is in the range of 10^ GeV. 
Only one large extra dimension is sufficient in our discussion. The smallness of the u^^^ is 
due to seesaw mechanism in the brane that uses intermediate seesaw scales (of the order of 
the string scale or so). The maximal z/^ — Uj. mixing needed for atmospheric oscillations owes 
its origin to the texture of the right handed neutrino matrix. The LSND data is explained 
in a natural way via small mixing between generations. Thus the new contribution of this 
paper is the identification of the lowest KK mode of the bulk neutrino as a viable sterile 
neutrino candidate and its embedding into a realistic four neutrino gauge theory framework. 

II. MODELS WITHOUT RIGHT HANDED NEUTRINOS 

Let us start by summarizing the simplest mechanism to produce neutrino masses using 
extra dimensions as discussed in Ref. [IU|. One assumes that all the Standard Model particles 



are localized on a brane embedded in the bulk of larger dimensions. The conservation of 
gauge flux then implies that, besides gravity, the only fields that could propagate in the 
extra dimensions are standard model gauge singlets. Their coupling with the brane fields is 
naturally suppressed by the volume factor -^^, where, from the observed strength of gravity 
the Plank scale and the string scale M* are related by 

M|, = M,"+V, (4) 

with Vn the volume of the extra space. We will assume that one of the radii of the compact 
extra dimensions, R, is larger. Now, by including the coupling of a bulk neutrino vb{x^, y) 
to the standard model lepton doublet L{x^, y = 0) we get the following terms which are 
responsible for the neutrino mass in this model 

S = K (fxLHi'BR^x, y = 0) + d'^xdyi'BLi.x, y)d^i'BR{,x, y) + h.c. (5) 

where H is the Higgs doublet, L is the lepton doublet and n is the suppressed Yukawa 
coupling 

By introducing the Fourier expansion of the bulk field we may write down the Dirac mass 
terms in fEI) as 







where u'^ embodies in a compact way the KK excitations along the fifth dimension and 
m = Kv is the mass term produced by the vacuum expectation value, v, of the Higgs field. 
The entry \/2Tn has to be interpreted as an (infinite) row vector of the form \/2m(l, 1, ■ ■ ■). 
The operator d^ represents the diagonal and infinite mass matrix of the KK modes. The 
eigenvalues of this operator are n/io, with /iq = 1/R. For ?ti -C /iq the mixing angle of the 
standard neutrino with the n-th bulk mode is [|l^] 

tan^„ ^ -; with ^ = V2mR. (8) 

n 

Then, the standard neutrino has a mass m and oscillations into bulk neutrino are present 
dominated by the first exited mode of mass /io- An explanation to the solar neutrino 
problem then follows if one assumes MSW mechanism with /xg — 10~^ eV {R ~ 0.2. mm). 
The question now is how one can extend the model to incorporate the atmospheric data and 
the LSND results. 

First thing to note is that in this scenario, the mixing between neutrinos of different 
generations is expected to be produced through the couplings with the bulk neutrino only. 
However, even if the bulk neutrino is not blind and couples with a different Yukawa coupling 
to each brane neutrino, the mass matrix will produce two massless neutrinos which decouple 
from the oscillation pattern. Therefore, there is no room for understanding the atmospheric 
and LSND data in this framework. 



Let us now see if this problem can be cured by including three bulk neutrinos. The most 
general Dirac mass in this case may be written, after a proper rotation of the bulk fields, as 



£. = vl- VMd ■ VBR{y = 0)+ I dy vbl ■ db ■ ^br + h.c. 



(9) 



with Ul = {i'e,i^fi,^T)L', i^B = {.^b^^bi^b)'i ^d = Diag{mi,m2,m^); and where the unitary 
matrix V depends on the texture of the Yukawa couplings to bulk neutrinos. The mass 
parameters in Mr, are small numbers of the order of eV or so and are given by the eigenvalues 
of the Yukawa matrix multiplied by v. Since the extra dimensional volume suppression is 
present, it is not unnatural to choose one of the masses to be in the range of 10~^ eV as in 
(^) to solve the solar neutrino puzzle. 

After rotating the standard sector by V , the last expression simplifies to the form 



C= Y. Wa^l^lRiv = 0)+ [dy u^BLdbi^BR + h.c. 

1 O O L " 



a=l,2,3 



(10) 



Therefore, in this picture the pattern of neutrino oscillations follows the same as before with 
a generational mixing given by V, which must provide atmospheric and LSND mixings. 
However a proper understanding of the neutrino puzzles would need three independent 
mass parameters. At first glance we would be tempted to believe that these masses could be 
those involved in (0). However, this is not the case. As we shall show below the heaviest 
eigenstate z/£ has a mass rha = ram{fio/2,ma}. Thus, if solar neutrino data is assumed to 
be solved by u^ — ub oscillation, this would require that /ii ~ 10^^ eV. As a result other 
mass differences in this picture become too small to solve the atmospheric and LSND data. 
If on the other hand, we assumed that it is the atmospheric neutrino puzzle which is solved 
by using the bulk neutrino, then we would choose fii ~ 0.06 eV; one could envision solving 
the solar neutrino puzzle by z/g — z/^ oscillation, with the relevant mass difference coming 
from small Yukawa couplings. Then we would be unable to explain the LSND data. 

Thus it appears that the only way out is to assume that the bulk neutrinos come from 
different branes with different sizes: /Xq — 10~^, and fi2,3 — 1 eV, with Am^^^ = /^2 — /^s- 
The challenge in such a scenario is to explain how the mass terms can be diagonalized 
simultaneously, a fact used to obtain Eq. (jTOl). 

To prove our previous statement, let us consider again the mass matrix (J^). In this 
compact notation it is simple to get the exact solution for the eigensystem even when the 
mass matrix is infinite. First the characteristic equation is given by 



{ml - dl) 



m„ — m + 






mt 



0, 



(11) 



where m„ is the mass eigenvalue and a sum on the last term is implicit. This expression 
translates into the same result that was obtained before in references HTsHTQI 



TTlr, 



nm^R cot{7imnR). 



The eigenstates are thus given symbolically by 



l^nL 



N„. 



l^L + 



ml - dl 



'BL 



(12) 



(13) 



where the norinahzation factor A^„ is easily computed to be 

A;2 1 



Nl = l + 2m^R^ J2 



Now we may express ul in terms of the massive modes just as 



\m J 



(14) 



OO -| 

^L = Y. W^nL- (15) 

n=0 ^^" 

Since Nn has a minimum for the lowest value of m„, the main component of vl is always 
the lightest mode z/ql. Notice from (|1^ that when mR <^ 1, the eigenvalues comes out to 
be rriQ = m and m„ = n/io for n > 0, and A^^ ~ mn/\/2m, thus, we recover the mixing angle 
(^). However, in the other limit, when mR ^ 1 the lowest eigenvalues are shifted down to 
{2n + l)fio/2. Then the main component of ul will have a mass of /io/2 and the mixing 
angle will be in general tan^ji = No/Nn- 

In consequence, irrespective of the hierarchy of the masses in Mo, a solution to the 
neutrino puzzles with bulk neutrinos needs localization of the three bulk neutrinos in different 
branes with extra assumptions about the Yukawa couplings. 

III. MODELS WITH RIGHT HANDED NEUTRINOS IN THE BRANE 

We now proceed to consider the second class of models where we include both the left 
and the right handed neutrinos in the same brane and a blind bulk neutrino which only 



couples with the right handed neutrinos due to either gauge symmetries |21[] or asymmetric 



boundary conditions in the bulk [23 



The simplest gauge model where this scenario is realized is the left-right symmetric model 
where the right handed symmetry is broken by the doublet Higgs bosons Xni^y 2, 1), where 
the number inside the parenthesis correspond to the quantum numbers under SU{2)l x 
SU{2)r X U{1)b-l- The relevant terms of the action for one generation are 

S = f (fx[nLxLVB{y = 0) + nRxRVB{y = 0) + hlcjyR] + j d^xdyvB^^d^VB + h.c. (16) 



As discussed in reference [^, by setting < x^r >= vr. and < xl >= O5 ^^^e following Dirac 



neutrino mixing matrix is obtained 

(^eL l^OBL 1^'bl) \ ^^R f ''f^ ] , (17) 

with pqbl being the zero mode and u'^ representing the exited modes as before. (Let us note 
parenthetically that in general a nonvanishing value for < x° > could be expected from 
the potential, so we should assume either that this vev is small so that its contribution to 
neutrino masses is negligible or that the bulk neutrino breaks explicitly the parity symmetry 
so the coupling with L is zero as advocated in Ref. [^.) Now, provided that kvji ^ /if ~ 



few MeV, a massless field which is predominantly the electron neutrino appears. Since 

6 



-^— , this constraint implies that M* must be as large as 10^ GeV or so. Oscillations 



into bulk neutrino will now result, again dominated by the lowest mass of the bulk modes pT 



implying that the largest radius of the extra dimensions should still be at the millimeter 
range. 

Let us now look at the complete picture for three generations. The more general mass 
terms in this class of models are 



C = Vl- MlL -Vl + Vl- Md-Vr + Vr- Mn -Vr 

+Vr ■ m ■ UBL{y = 0) + I dy Ubl ■ <9; 



k ■ ^BR + h.c. 



(18) 



In last equation M.ll represent the Majorana mass terms of the left handed neutrinos. For 
simplicity we will assume those terms to be zero. M^ represents the Majorana matrix for 
the right handed neutrinos which in the left-right models arises from the vev oi a. B — L = 2 



SU{2)r triplet and m 



m 



1, 1, 1) are the universal couplings of the bulk neutrinos to the 



right handed neutrinos. The mass matrix now has the profile 



V 



Md 








\ 


Ml Mn 


m 


v/2m 





mt 











V2mT 








^5 








^5 


oy 



(19) 



in the basis {vli1-'Ri1-'obli1-''bli^'bii)- Notice that this is really an infinite matrix with the 
entry \/2m. being the universal coupling of ur to the exited bulk modes VuBL-, then, this is an 
infinite matrix itself given by -\/2m x (1, 1, ■ ■ ■) . We have used the two component notation 
and used subscripts L, R as labels for corresponding states rather than helicity projection 
operators. Notice that fosR is massless and is therefore not included in the above mass 
matrix. As a result, this matrix has odd number of rows and columns and has one zero 
eigenvalue. The corresponding state is identified due to its fiavor content to be the close to 
the z/g state. The elements of Mat are expected to be large, < M=k, while those in Md as 
well as m are of the order of some MeV. This via the usual seesaw mechanism leads to light 
eigenstates. 

We can write down the general mixing matrix for the "four neutrino" states as: 



\^bJ 



u ■ 






(20) 



where the neutrinos in the right hand side are the mass eigenmodes, and v^ represent the 
mass eigenstates of the KK tower. The active neutrinos will contain in general suppressed 
contributions from the complete KK tower, then, only the lowest mode is expected to con- 
tribute substancially and we will at the end identify this mode with v^. We will show that in 
this picture, \Uei\ ~ 1 and \Uei\ > |f^e2,3|; ^i being the lightest element of the spectrum with 
zero mass. If the masses of z/2^3 are in the eV range with a mass difference square of the the or- 
der of Aml^^, we can solve the atmospheric neutrino puzzle if | ?7^2 1 / 1 f^/^s | ~ |f^r2|/|f^r3| ~ 1, 



with the decouphng of the bulk modes from this sector; thus we would need |f/^,T4| <^ 1 and 
\UBi\ ~ 1- Finally to explain the LSND data, we need \U^i\ ~ \Ue2\ 3> |t/^e3,/i4|- We will 
show below that all these features emerge only from the simple choice of the seesaw matrix 
Mjv as is commonly done in most model building [Q . 



Let us see how this scenario is realized in our picture. First note that the mass matrix 
in ( p!9D has a zero mass left handed eigen mode, un, which is a linear combination of ul and 
uoBL given by the solution to the matrix equation 

[MJ,m)-\iyiL>=0. (21) 

where ( M^j m ) is a 4 x 3 matrix. If we want this lightest eigenstate to be predominantly the 
electron neutrino to maintain the observed universality in charged current weak interactions, 
we must demand that in the expression for un, 

^IL = -T7 WeL + Siy^L + ^l^rL " et'Oi] (22) 

we must have (5, A, e ^ 1. If we make the simplest choice for M^ = Diag(mi, m2, m^), with 
the natural hierarchy mi < m2 < rris, we get 6 = mi/m2 > A = mi/m^ and e = rrii/m. We 
then may set 1 > e > 5 3> A to get sin^ 29soi — e^, and S is going to produce the mixing for 
the LSND data. 

After the extraction of the zero mass term ( P^D we may procceed to decouple the right 
handed fields which get large masses via Mj^i. The mass matrix for light fields then reduces 
to 

MlM^^mI V2MiM^^m 
\/2mt M^^M^ 2mtM;v^m ^5 | , (23) 

^5 

in the basis {vqbl, ^(iL, ^tL, ^'bl^ ^'br)i with Ml of the form 

m m m \ 



Ml= \ m2 

mg y 



(24) 



where we have neglected those contributions that are proportional to the small parameters 
e, 5 and A to show how the decoupling of z/^,,- may occurs. It is worth noting that this 
form of Ml is produced by the conservation of L in the Standard Model which suggest that 
Mo should be diagonal, thus giving Ml as above. In this picture, the L violation terms 
come from high energy, and all the texture in the mixing matrix is produced through M^- 
Also, notice that the mass matrix ( ]23|) is really of a generical form. It may be produced not 
only by see-saw as we are assuming here but as well as through some other mechanism as 
radiative corrections. 

Now, the main ingredient of the decoupling is based on the fact that the bulk neu- 
trino really couples to a certain linear combination of the brane neutrinos, given by 
MLM^^iyoB, i^n, ^t)l^ ^^d decouples from the orthogonal combinations. Moreover, the first 
row and column of MlM^^^M]^ correspond to the same combination. Therefore by rotating 
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the brane sector to diagonalize MlMj^^M^, we may expect that the decouphng becomes 
exphcit. Lets see how this argument works. First notice that if we make the rotation 



1^3 



K 






(25) 



then the mass matrix (p3|) transforms into 






KMlM 

V2 



y/2 iTMiM^^m 
2 m1'Mr/m 







d. 



N 



(26) 



Let us choose K such that K M^M^^ MJ^K^ becomes diagonal. All the mixing with the bulk 
modes i^'bl r ^^^ then given by the form of the column vector KMlM^^ui. It is here that 
we expect to see the decoupling. Let us take an explicit form for M]s[ to make it clear: 



M, 



N 




(27) 



and to simplify matters let us assume that m = 10, 1712 = 50, mS = 100 and a = 10^, 
6 = 5 X 10^ and c = 10® all in MeV. These values are selected merely for illustrative purpose. 
Other values may also achieve the desired decoupling. The correct entries in the matrix 
depend on the high scale sector of the model, which we leave unspecified for our purpose. 
Now we proceed with the numerical computation with above choice of numbers and get 



K 



0.9998 
-1.4- 10 
-1.4- 10 



2 X 10"^ -2 X 10^2 
-2 0.7103 -0.7038 
-2 -0.7039 -0.7002 



(28) 



while the eigenvalues of MlM^^mI are given by (.25 MeV, -1.0028 eV, 1.0008 eV ). Thus 
we get a maximal mixing in the right sector. On the other hand, for the mixing with the 
bulk modes we get KM^M^^m = (.25, 1.4 ■ 10-^ -1.4 ■ 10"®) MeV. Thus, 1/2,3 decouples 
while v'qbl mixes maximally to the bulk modes. 

Next, as the entry m'^'Mj^^in in (|23| ) also is the first element in the mixing MiM^^m^ 
after extracting the decoupling modes, the effective mass matrix reduces to the profile 



a y/2a 
(95 



/ a V2a V2a 

^/2a 2a /io 2a 



2a 2/io 

2/io 

V I I ; ; ; 






/UO 





'2a 


2a 















(29) 



where a = m^M^^m. It is simple to check this fact in the case above since a ~ .25 
MeV too. It is worth noting that now all the fields involved are bulk modes. The rows 



and columns in the left hand side of (^) are labeled by {i^oL,J^BLy^BR)- ^^ the right 
hand side we have explicitly written it in its true infinite mass matrix form labeled by 

i^OL, ^IBL, ^IBR, ^2BL, ^2BR, ' ' ■)• 

We may worry about the high value of a since it could mean that i'obl has a component 
with a large mass; however, as a result of the maximal mixing with the KK modes, that 
contribution will be suppressed as well as the other heavy elements of the tower are. The 
mechanism works exactly as in the one generational case of this class of models |^ where 
the ul which gets mixed with the bulk modes has a mass term in the MeV range [see Eq. 
(P^)]. There, the heavy term is absorbed by the tower and as a result the masses of the KK 
modes are shifted down, while their contributions to z/^ becomes suppressed as 1/n. This 
suppresses the heavy field contribution by the number of levels below such a mass i.e. ~ aR 
MeV, which is a large number. The lowest mode then becomes the main component of uql 
that appears in Eq. (21); We identify this mode as the sterile neutrino, 1/4, used in solving 
the solar neutrino puzzle. 

Let us demonstrate this by diagonalizing the mass matrix (^). We take advantage of 
this compact notation to write the characteristic equation as 



(A 



' dl) 



a 



A + 



2aA2 



A2 



dl 







(30) 



which now reduce to 7rQ;-Rcot(7rA-R) = 1 and may be solved exactly to get the mass eigen- 
values 



A. 



arctan(7rQ;i?) n 



(31) 



where now n E Z, since also the right handed bulk fields are involved. Clearly, when aR ^ 1 
(as is the case here) the masses are just A„ = {2n + l)/io/2, which are now shifted down by 
\lXo. Also we note that the only masses of order of a correspond to n ~ ±ai? which were 
already in the tower; thus our massive term has been totally absorbed into the KK tower. 
Strictly speaking, there is an extra mode in the tower with infinite mass. 
The corresponding eigenvectors are given by (A are the eigenvalues) 



KA) 



1 



with the normalization factor: 



l^OBL + 



x/2A2 



i^m^ + 



v^A(95 



A2 _ 92 '^BL ^ ^2 _ Ql 



'BR 



(32) 



^ 



1 + {TxaRf ^ (2n + 1 



TT 



(33) 



Notice that in our limit i.e. aR ^ 1, the masses are degenerate, and occur in pairs with 
masses ±An. We may then recombine those states to get two Majorana neutrinos that form 
a Dirac neutrino of mass A„ of the form 



VnL 



V2 

^(An^ 



VQBL + 






di 



V 



BL 



1 



^nR 



2Xnd^ 



ViK) X'r 



a2 



'BR 



(34) 
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(36) 



Now, we may rotate the system backwards to write down uqbl in terms of the left handed 
mass eigenstates to get 

Thus, as rj grows quadratically in A„, the main element of z/qbl is the lowest mode z/ql that 
we then identify as the sterile neutrino. Moreover, from (|33D we see that the suppression 
over the other modes goes as l/(2n + 1), then the heavier modes get highly suppressed, since 
n is also the number of eigenlevels below. 

Finally, we may put all the ingredients together and calculate the effective mixing pattern 
at low energies. We have proceeded using a numerical analysis from the beginning to estimate 
masses and mixings directly from the mass matrix in (|19|), without the approximations we 
made in the former analysis to get more accurate values. The input values for our parameters 
are the same as before with mi = 1 MeV and we have taken /^o = 5 x 10^^ eV to fit properly 
Am'j.^i , we then get 

/ z/i(massless) \ 

( i^e\ _ ( 0.9948 0.0212 -0.0072 0.0632 \ z/2(-l-0014 eV) 
\i^^J ~ \ 0.0199 -0.7037 0.7102 0.0013 ) 2/3(1. 003 eV) 

V z/4(2.5151 ■ 10-3 eV) J 

where the masses are indicated between brackets. From these results we get 

Am'i^i ~ 6.326 x 10-*^ eV^, sm'^20soi ~ 1.58 x lO'^ 

Aml^^ ~ 3.23 X lO'^ eV^, sin"^ 29atn, ~ .9991 (37) 

Aml^ ~ 1.0028 e\/2^ sin^ 2^e/. -1.18x10^3 

These values are to be compared with the experimentally preferred values ([l|), (H) and @. 
It clearly shows that the choice of the bulk neutrino as the sterile neutrino works quite well 
to explain observations. Although for the sake of explicit demonstration we have made a 
specific choice of parameters, there is a range of parameters in the mass matrix that will do 
equally well. 

IV. CONCLUSIONS 

In conclusion, our analysis seems to demonstrate that one can use the bulk neutrino as a 
sterile neutrino needed to understand all neutrino oscillation observations. Despite the fact 
that it has an infinite set of KK excitations, only the lowest mode seems to play an impor- 
tant role and indeed gets identified as the sterile neutrino. It is worth emphasizing again the 
conceptual advantage of this approach which is that the ultralightness of the sterile neutrino 
has a geometric origin rather than from extra symmetries, as in the four dimensional mod- 
els. What is also gratifying is that the apparent connection between submillimeter gravity 
experiments and the neutrino puzzles seems to remain in the realistic implementations of 
this idea. 
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